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MAPPING THE GALACTIC HALO WITH BLUE HORIZONTAL BRANCH STARS FROM THE 2DF QUASAR 

REDSHIFT SURVEY 
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ABSTRACT 

We use 666 blue horizontal branch (BHB) stars from the 2Qz redshift survey to map the Galactic halo in four 
dimensions (position, distance and velocity). We find that the halo extends to at least 100 kpc in Galactocentric 
distance, and obeys a single power-law density profile of index ~ —2.5 in two different directions separated 
by about 150° on the sky. This suggests that the halo is spherical. Our map shows no large kinematically 
coherent structures (streams, clouds or plumes) and appears homogeneous. However, we find that at least 20% 
of the stars in the halo reside in substructures and that these substructures are dynamically young. The velocity 
dispersion profile of the halo appears to increase towards large radii while the stellar velocity distribution is non 
Gaussian beyond 60 kpc. We argue that the outer halo consists of a multitude of low luminosity overlapping 
tidal streams from recently accreted objects. 
Subject headings: Galaxy: formation — Galaxy: halo — stars: horizontal-branch 
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1. INTRODUCTION 

The motion of old stars preserves the fossil record of the 
earliest phases of formation in the Milky Way, as was realized 
in the two seminal papers by Eggen, Lynden-Bell & Sandage 
(1962) and Searle & Zinn> (119781) . The picture that has since 
emerged is one where the Gala xy has been built via a series 
of acc retions and mergers (e.g.. iFreeman & Bland-Hawthornl 
2002), in agreement with hierarchical structure fo rmation sce- 
narios (Johnston et al. 2008; Cooper et al. 2009). Consistent 
with these theoretical expectations, the halo is thickly popu- 
lated by debris from past and on-going accretion events (e.g., 
Ibata, Gilmore & Irwin 1994; Belokurov et al. 2006). Promi- 
nent debri s features have also been observed in the halo s 
of M 3 1 dlbata et all |200H 1200 7; McConnachie etail l2009h . 
M33 dlbata et al. 2007; McConnachie et al. 2009) and in other 
nearby spirals (Martinez -Delgado et alJl2010h . 

At the same time, the inner halo (at galactocentric dis- 
tance Rgc < 20 kpc) of the Milky Way contains an old, 
dynamically smooth and metal-poor component, which was 
proba bly formed by rapid early merging or vi olent relax- 
ation (iCarollo et alJl2007i 120091 feell et alJl2008h . While the 
buildup of most of the inner halo is believed to take place 
rapidly, the outer regions of the Milky Way should grow more 
slowly via the disruption of dwarf galaxies: the outer halo is 
therefore expected to be quite inhomogeneous and dominated 
by large streams, clouds and plumes from infalling objects 
(Johnston et al. 2008; Cooper et al. 2009). The main question 
we need to address is the relative role of mergers and accre- 
tion vs. in situ formation, i.e., whether the observed struc- 
tures represent the main mechanism by which the Galaxy was 
formed or are only a comparatively minor contribution over 
an ancient stellar component. 

In order to clarify these issues we explore the space dis- 
tribution and kinematics of the outer halo using Blue Hori- 
zontal Branch (BHB) stars. These objec ts have often been 
used as tracers of the Milky Way halo (e.g., Yanny et al. 2000; 
Kue et al. 2008; Bro wn et all llOlO). BHB stars are compara- 
tively bright and reliable standard candles, are more common 
than other commonly used probes (such as Carbon stars or RR 
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Lyrae) and are representative of the old and metal poor stellar 
population that constitutes the Galactic halo. These stars can 
be easily identified spectroscopically, even at comparatively 
moderate resolution, but are sufficiently rare that extensive ra- 
dial velocity surveys are needed to obtain a significant sample 
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(lYannv et al.ll2000[|Xue et aT]l2lMlBrown et al.ll20Toh . 

Here we exploit a large suite of archiv al spectroscopy fr om 
the 2dF Quasar Redshift Survey (2Qz - iCroom et a l. 2004) to 
identify a sample of 666 BHB stars in the halo out to 100 kpc 
and explore its structure and kinematics. We use this dataset 
to set limits to the size of the Milky Way halo, its shape, pres- 
ence of streams, degree of substructure and the velocity dis- 
persion profile out to large radii. The next section describes 
the 2Qz survey and how we identified BHB stars. Section 3 
discusses the 4D map of the Galactic halo we produce from 
these data. Section 4 examines the question of substructure 
in the halo and the final section analyzes the halo kinematics 
and discusses the results in the context of galaxy formation 
models, especially as apply to the Milky Way. 

2. THE 2QZ DATA: IDENTIFICATION OF BHB STARS 

The data for this project come from archival observations 
of the 2Qz survey. The 2Qz obtained spectra for ~ 50, 000 
A-colored point sources with 15 < bj < 20.9 selected from 
Supercosmos survey photometry (Hambly et al. 2001J). The 
targets lie in two 75° x 5° strips on the sky, one on the celestial 
equator between 09h 30m < RA < 14h 30m and -2.5° < 
5 < +2.5° and the other centred in the Southern Galactic 
Cap between 22h 00m < RA < 03h 00m and -32.5° < 
5 < —27.5°. We refer to these as the Northern and Southern 
samples, respectively. 

The targets were selected on the basis of their u, bj and rp 
colors (photographic phot ometry from the SuperCosmos sur- 
vey - lHamblv et alj|200ll) t o lie in the param eter space cov- 
ered by known quasars (see[Croom et al. 2004 for further de- 
tails): targets were selected to have bj < 20.9 and —2.5 < 
u — bj < 1.5 and— 1.5 < bj — rp < 2.5, excluding a box with 
colors 0.4 < u — bj < 1.5 and < bj — rp < 2.5 which con- 
tains main sequence stars. Spectroscopy for these objects was 
obtained during 1998-2003 with the Anglo-Australian 3.9m 
Telescope in Coonabarabran, NSW, Austr alia, using the 2 de - 
gree field (2dF) multi-object spectrograph dLewis et al.l2002h . 
The spectra cover the entire optical range (3500 < A < 7500 
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FIG. 1 . — The color distribution of BHB stars identified in the 2Qz survey. 
The original Supercosmos colors have been transformed to the SDSS using 
stars in common between the two surveys. 

A) at a resolution of about 2000, with exposure times of at 
least one hour per target, although a fraction of the objects 
were exposed for considerably longer, if they happened to lie 
in a region where the 2dF tiles used by the survey overlapped 
(to insure greater completeness) or the field was reobserved 
(because of lower than expected signal). 

As with all color-selected surveys, 2Qz suffers from signif- 
icant contamination from QSO-colored stellar objects, such 
as white dwarfs, disk A stars, blue stragglers and BHB stars, 
although this is minimized by observing at high galactic lat- 
itude. The data were reduced and redshifted via a semi- 
automated technique and whenever a stellar redshift (z = 0) 
was obtained, the object was discarded from further analysis 
but placed in the database. We proceeded to retrieve all stellar 
spectra from the database and classify them on the basis of the 
equivalent widths of the H 7 and H^ Balmer line s, to id entify 
a samp le of 666 bona fide BHB stars. Following lYanny et al] 
(2000); IXueetalJ d2008l) : iBrown et alj ( 120101) we first mea- 
sured radial velocities for all star by c ross correlating with 
synthetic templates from the library of Mu nari et alj (J2005): 
the templates had temperatures and gravities typical of A-stars 
and field horizontal branch stars dBeers e t al. 2001). We then 
fit Gaussian curves to the H7 and H<5 lines and measured the 
width of the Gaussian fit at 20% of the normalized contin- 
uum level: this indicator Do. 2 has been shown to be a good 
discriminator betw een BHB stars, blue stragglers and other 
contaminants (lPierlll983t lYanny et al. 2000). We only used 
spectra that we deemed to be of sufficient quality to allow a 
secure classification. In order to be included in our sample ra- 
dial velocities had to be determined to within 50 km s _1 and 
the H7 and H<5 widths had to have errors of less than 20%. In 
order to classify stars as bona fide BHB stars we require that 
the mean £>o.2, from both lines, lie between 17 and 31 A (as 
in lPierlll983tlYanny et alJl2000 . leaving a total of 666 BHB 
stars in our sample. 

Figure 1 shows the distribution of the stars in the u — g vs. 
g — z plane, where we transformed our Supercosmos u — bj 
and bj — rp colors to the Sloan system by using stars in com- 
mon in the equatorial region shared by 2Qz and the SDSS. 
This is at least qualitatively similar to the color distribution 
of sta rs in previous work (e.g., compare Fig. 1 in Brown et al. 
2010) and suggests that our sample is comparable to those 
used in previous studies, in terms of selection criteria and 
degree of contamination from blue stragglers and other A- 
colored objects. 

As our targets span a relatively narrow range in colors (as 



selected by the 2Qz survey) we assumed an absolute mag- 
nitude of M b = .7 ± 0.2 which is typical for BHB stars 
dLavden et al.l [T996). We finally used these distances, the 
known sky positions and the measured radial velocities to 
place all our stars on a cylindrical coordinate system at rest 
with respect t o the centre of the Galax y, assuming Solar 
positions as in iDehnen & Bin nev (1998). This yields a 4- 
dimensional (position, distance and radial velocity) map of 
the galactic halo in two widely separated (150°) lines of sight. 

3. THE 4D STRUCTURE OF THE MILKY WAY HALO 

In Figure 2 we plot the 4D map of the Galactic halo we pro- 
duced, projected along the three most relevant dimensions. It 
is clear from this figure that the halo of the Milky Way extends 
to at least 100 kpc from the Galactic centre, and likely well be- 
yond (the edge of the map is set by the magnitude limit of 2Qz 
data), in both directions we survey. This is considerably larger 
than previously believed and comparable to the large metal- 
poor halo observed in the Andromeda galax y dChapman et alj 
l2006t lKalkailFaT]|2006l iKoch et alj|200g) : it would include 
several of the dwarf galaxy satellites (including the Magel- 
lanic clouds) within the Galaxy's stellar halo. As a matter of 
fact the Sextans dwarf is visible in Fig. 2 at x ~ —20 kpc and 
y ~ +50 kpc. Such large metal-poor halos may be ubiqui- 
tous (e.g., in NGC 3379 - lHarriset all 120071) and may be a 
common byproduct of early galaxy formation. 

In previous work, the SDSS has identified BHB stars out to 
60 kpc from the Galactic centre (Xue et al. 2008), while the 
Hypervelocity Stars Survey (IBrown et al . 2010) found a BHB 
star sample to a distance of 75 kpc. The S paghetti survey 
dMorrison et all 120001: IStarkenburg et alj|2009l) has observed 
halo red giants to a distance of 100 kpc, and star counts in 
the COSM OS field identify a halo component to a distance 
of 80 kpc (Robi net al.ll2007l) . Our study reaches deeper than 
nearly all these and covers a larger field of view than all except 
the SDSS (it is comparable to the Hypervelocity Stars Survey 
coverage). However, we sample the BHB stars more densely 
as all potential targets have been targeted by the 2Qz survey, 
although of course we are not able to classify all stars. Since 
we cover nearly diametrically opposite areas on the sky, we 
argue that the detection of the stellar halo in our data is not due 
to possible diffuse structures on the sky that accidentally lie in 
our line of sight (as for pencil-beam studies such as COSMOS 
or the Spaghetti survey) and the Milky Way halo truly extends 
to large radii. 

Figure 3 shows the radial density profile of the halo along 
both directions we survey. In both cases we obtain a good fit 
to a single power law of index i?~ 25 ±02 . This is somewha t 
shallower than the ~ R~ 3 found by iMorrison et al.l d2000l) 
and predicted by theory, but is in good ag reement w it h pre- 
vious measurements using BHB stars by IXue et ail d2008l) 
and Brown et al. (2010). We are of course incomplete in that 
we cannot detect and identify all BHB stars in 2Qz. This 
incompleteness is a complex function of our ability to reli- 
ably classify stars as a function of spectroscopic signal-to- 
noise. Naively, we would preferentially miss the most dis- 
tant objects, that would tend to make the radial profile steeper 
than it actually is, while contamination from blue stragglers 
would tend to make the profile flatter. T he similarity betwee n 
our profiles and thos e derived by (e.g.) lYanny etal J (2000); 
IXue et all d2008l) and lBrown et alj d2010l) argues that our con- 
tamination fraction (from blue stragglers) and completeness 
are not too different from previous samples. 

In any case, this should not affect our differential measure- 
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FIG. 2. — The 4D map of the Galactic halo, with stars projected on the Galactic plane (along a line connecting the Sun to the Galactic centre) on the x-axis and 
vs. their height above the plane on the z-axis. The stars are color-coded by their Galactocentric radial velocity, as indicated in the figure legend. 
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FIG. 3. — Radial density profiles for BHB stars in the halo and best power- 
law fits to the Northern and Southern samples separately. 

ment of the radial density profile for the two individual sight- 
lines we survey. The similarity in the profile slope then argues 
that the stellar halo of the Galaxy is spherical (Majewski et al. 
120031: Smith, Wyn Evans & Ahn 2009), althoughmore sight- 
lines would be helpful to obtain a more precise estimate for 
the halo's axial ratio. 

In a recen t surv ey of RR Lyrae in SPSS stripe 82, 
IWatkins et a l. (2009) and lSesar et all (1201 Oh find that the halo 
radial profile becomes very steep (i?~ 6 ) at galactocentric 
radii larger than 40 kpc and the the outer halo seems to be 
dominated by a large structure known as the Virgo-Aquila 
cloud, which is either an infalling dwarf galaxy or a tidal 
str eam. However this is not observed in our data, or those 
of iBrown et al.l d2010l) . One possibility is that the RR Lyrae 
distribution is more concentrated towards the Galactic centre, 
as RR Lyrae tend to be more metal rich than BHB stars and 
the halo appears to have an abundance gradient (Carollo et al. 
2007, 2009). This would create an apparent break in the radial 
profile, reflecting the lower overall metallicity of the sample 
at large distances. 



4. SUBSTRUCTURE IN THE HALO 

Inspection of Fig. 2 also shows that there are no obvious 
large kinematically coherent features, such as streams, plumes 
or shells, in our 4D map of the Galactic halo. We identify 
three possible structures: one is the Sextans dwarf, as pre- 
viously mentioned. A plume of infalling stars is present in 
the top right-hand corner of Fig. 2 at x ~ +40 kpc and 
y ~ +80 kpc. This may be related to the Virgo-Aquila struc- 
ture. Finally, a small clump of objects is observed in the di- 
rection of Piscis Austrinus at x ~ +40 kpc and y ~ —50 kpc 
and may be an undiscovered dwarf. Nevertheless, we do not 
appear to observe the large stream s expected in simulations 
dJohnston et alj|2008t ICooper et aT]|2009l) and the outer halo 
appears to be more homogeneous than predicted. 

We can quantify the presence of streams or oth- 
erwise using the Great Circle Stream method of 
iLynden-Bell & Lynden-Belil (119951) in Figure 4. We 
compute the excess radial energy of each star compared to 
a smoothed Galactic potential and compare the results to a 
randomized sample, where we scramble the velocities, but 
not the positions, of stars, 1000 times. Streams would show 
up in this figure as long lines of stars. It is easy to see that 
our sample contains no significant stellar streams. The large 
metal poor outer halo appears to represen t an extension of the 
smoot h metal poor in ner halo structure (ICarollo et al.ll2007l 
l2009l:lBeiretaLll2008h . 

Howeve r, lack of large streams does not mean lack of sub- 
structure. iBell et al.1 (12008) find that about 40% of stars in 
the inner halo are substructured, even if the only prominent 
stream in their data is t he we ll known Sagittarius stream, 
while IStarkenburg et alj d2009l) find that around 25% of the 
stars in the Spaghetti survey are more clustered than expected 
from a random distribution. We ap ply the 4-distance method 
used by Starkenbur g et alj d2009l) to our data in Figure 5. 
This is essentially a version of the correlation function for 
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FIG. 4. — The Great Circle Stream method applied to stars in our sample. 
We plot the excess radial energy of each stars with respect to a smoothed 
Galactic potential vs. its radial distance from the centre of the Galaxy. The 
left hand panel shows the actual data, where the right hand panel is a random- 
ization. No streams are observed in these figures as long lines of stars having 
the same excess radial energy. 
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FIG. 5. — The excess fraction of stellar pairs over random (see text) in 4- 
dimensional distance for our sample of BHB stars. 

4-dimensional data, computing the excess number of stellar 
pairs in position and velocity compared to a random sample 
(realized by scrambling the velocities but not the positions of 
the stars in the data) as a function of distance in 4-dimensional 
space. The angular, spatial and velocity separations are scaled 
by the range in these quantities covered by the data. 

At a 5cr level we find that at least 20% of our stars are more 
pa ired than random. This agrees, broadly, with the estimates 
by Bell et al. (2008) for the inner halo and Starkenburg et al. 
(2009) for the outer halo. In addition we detect a decrease in 
the correlation strength at small 4-distance. This is expected if 
the outer halo is dynamically young and suggests the presence 
of numerous streamlets and a complex structure, too weak to 
be resolved by our data. We present further evidence to this 
effect from an analysis of halo kinematics. 

5. KINEMATICS OF THE HALO 

The kinematics of metal-poor stars in the halo yield 
information on th e earliest phases of galaxy formation 
(Egge net al.lll962l) and the shape of the Galactic potential. 
Until recently, known samples of halo s tars were small, espe- 
cially at large galactocentric distances. Battaglia et al. (2005) 
used a heterogeneous sample of red giants, BHB stars, glob- 
ular clusters and dwarf galaxies to trace the velocity disper- 
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FIG. 6. — Top panel: Run of radial velocities of BHB stars as a function 
of Galactocentric distance. Bottom panel: The radial velocity dispersion pro- 
file of the Galactic halo in both Northern and Southern samples. Each bin 
contains equal numbers of stars and the velocity dispersion and its error are 
calculated using a maximum likelihood method. The error bars on points for 
the x-axis represent the size of the bin used. 



sion profile of the Galaxy out to ~ 100 kpc , finding a mildl 
declining profile. Using BHB stars in SDSS Xue et al.1 (1200: 
found a flat or mildly declining profile out to 60 k pc and there- 
fore in ferred a large mass of the Milky Way. iBrown et alJ 
d2010l) also derive a mildly declining profile out to 75 kpc us- 
ing BHB stars in the Hypervelocity Star Survey. In all these 
cases, the fits depends strongly on the accuracy of the last 
(most distant) data points, which generally contain fewer ob- 
jects. 

We use our data to derive the velocity dispersion profile in 
both Northern and Southern samples separately. We use bins 
containing equal numbers of stars and calculate the velocity 
dispersi on and its error fol lowing a maximum likelihood ap- 
proach (IWalker et al. 2006). Figure 6(a,b) plots our results. 
While we are in reasonable agreement with previous work 
over the range where we overlap, we find evidence of a rising 
velocity dispersion at large radii, in both fiel ds w e study. This 
is unli ke the results of Battagli a et alJ d2005h and lBrown et alJ 
d2010l) . although we reach farther into the distant halo than 
they do. In the former case, the heterogeneous sample and 
its relatively small size may cause part of the difference, as it 
is kno wn that differ ent tracers have different kine matics (e.g., 
Kinmanet al. 2 0071) . In the latter case, we and IBrown et alJ 
(2010) use the same tracers , but while our star s are concen- 
trated in two narrow strips, IBrown et alJ (120101) uses a wide 
region selected from the SDSS. We have more stars (by a fac- 
tor of about 3) than they do at large (R > 50 kpc) distances. 
One possibility is that we sample this regime more densely, 
especially if the halo is as inhomogeneous as theory suggests 
and as the evidence from Fig. 5 also argues. 

Is it possible that our result may be affected by contamina- 
tion from blue straggler stars. Since these objects are 2 to 3 
magnitudes fainter than BHB stars, blue stragglers will con- 
taminate the bins at Rgc > 40 kpc with objects that truly lie 
at distances of 10 to 30 kpc. The selection procedures we de- 
scribe above should exclude most contaminants but even if we 
are as effective as previous studies, ~ 20% of our sample may 
consist of blue stragglers. In order to assess the effects of con- 
tamination, we have used a model distribution function for the 
velocity dispersion profile of our Galaxy from van Hese, Baes 
& Dejonghe (2009), with Galactic mass of 1.1 xl0 13 Mfpan d 
halo scale radius of 40.5 kpc from lDehnen & Binnevl dl998). 
We sampled 100 stars at distances of 10 to 100 kpc with a 
Gaussian random deviate having a(r) from Ivan Hese et alJ 
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FIG. 7. — The effect of BHB star contamination on the velocity dispersion 
profile of the Galaxy. The top panel is for a distribution with no velocity 
anisotropy while the bottom panel assumes that all velocities are tangential 
at large radii. See text for details. The degree of contamination assumed is 
indicated in the legend. 



(2009) and zero mean velocity. We adopted two cases: one 
with no velocity anisotropy /3q = j3i n { = and one with fully 
tangential radial velocities at infinity (/3o = 0; /3i n f = 1). For 
each of these we adopted a variable contamination, between 
10% and 40%, from blue straggler stars, which we assumed 
to be 2.5 mag. fainter than BHB stars. We replaced this frac- 
tion of stars in each bin with Roc > 40 kpc, with blue strag- 
glers at the appropriate distance, with velocitie s sampled from 
a Gau ssian random deviate with cr(r) from Ivan Hese et al.1 
(2009), at the appropriate r for a contaminating blue straggler. 
The results of this exercise are shown in Figure 7. Although 
the blue stragglers increase the measured velocity dispersion 
at large radii (compared with the theoretical profile), they do 
not produce a flat or rising velocity dispersion profile (as in 
the bottom panel of Fig. 6) even at 40% contamination, and 
even for the relatively flat <j(r) profile produced by the fully 
isotropic model. 

In Figure 8 we plot the histograms of radial velocities for all 
stars in 6 bins, containing equal numbers of objects and cov- 
ering a range of distances. We see that while the distributions 
are acceptably Gaussian within the inner 60 kpc, we observe 
both an increase in dispersion and a loss of Gaussianity in the 
two outer bins. This is consistent with what we observe in 



Fig. 5 and 6, where we find indications of an increasing ve- 
locity dispersion and a dynamically young and substructured 
halo. If the outer halo consists of a myriad of streamlets, ac- 
creted more or less recently from the disruption of satellites 
or low luminosity galaxies, we would expect evidence of dy- 
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FIG. 8. — Histograms of the velocity distribution for 6 bins containing equal 
numbers of stars (the distance ranges sampled are indicated in figure legends). 
The two more distant bins appear to have a much less Gaussian distribution 
than the four bins at R < 60 kpc. 

namical youth and tidal heating, as we observe in Fig. 4 and 
5. 

A possible alternative explanatio n is that the dark halo is 
very massive. Gned in et al.1 (120 10J) has recently argued for a 
very massive halo based on the flatness of the velocity dis- 
persion profile. However, in this case the mass to light ratio 
needed would exceed several hundred. Tidal heating may also 
produce an increasing velocity dispersion, as is observed in 
some dwarf galaxies (Munoz, Majewski & Johnston 2006). 
These explanations, while possible, would not probably pro- 
duce the non Gaussian velocity dispersion observed at large 
radii. 

The combined evidence from the space distribution and 
kinematics of BHB stars points to a very large and metal poor 
halo, whose outer regions appear to have been accreted rel- 
atively recently from low luminosity satellites, analogous to 
the copious tidal debris observed in M31. Although large 
streams are not observed the data appear to be in compara- 
tively good agreement with theoretical predictions, although 
it is possible that minor mergers are more important than ex- 
pected. More observations as well as models to truly 'ob- 
serve' the simulated halos in the same fashion as the data and 
directly compare theory and reality will be needed to refine 
our understanding of the formation of the Milky Way. 
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